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Abstract 

Velocity and density structure factors are measured over a hydrodynamic range of scales in a 
horizontal quasi-2d fluidized granular experiment, with packing fractions (p G [10%, 40%]. The 
fiuidization is realized by vertically vibrating a rough plate, on top of which particles perform 
a Brownian-like horizontal motion in addition to inelastic collisions. On one hand, the density 
structure factor is equal to that of elastic hard spheres, except in the limit of large length-scales, 
as it occurs in the presence of an effective interaction. On the other hand, the velocity field shows 
a more complex structure which is a genuine expression of a non-equilibrium steady state and 
which can be compared to a recent fluctuating hydrodynamic theory with non-equilibrium noise. 
The temporal decay of velocity modes autocorrelations is compatible with linear hydrodynamic 
equations with rates dictated by viscous momentum diffusion, corrected by a typical interaction 
time with the thermostat. Equal-time velocity structure factors display a peculiar shape with a 
plateau at large length-scales and another one at small scales, marking two different temperatures: 
the "bath" temperature Tf,, depending on shaking parameters, and the "granular" temperature 
Tg < Th, which is affected by collisions. The two ranges of scales are separated by a correlation 
length which grows with (p, after proper rescaling with the mean free path. 

PACS numbers: 45.70.-n, 05.40.-a, 61.43.-j 
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I. INTRODUCTION 



Understanding granular systems is strategically relevant for the optimization of many 
industrial processes in pharmaceutical, food, cosmetic, chemical, petroleum, polymer and 
ceramic industriesii^. It is also essential for modeling geophysical phenomena such as sedi- 
ment fiuidization in landslides^"- and soil mechanics^. Finally, it is a challenge in the field of 
statistical physics of out-of-equilibrium complex materials^"—. Experimental granular ma- 
terials are usually limited to a few thousands particles, sometimes even less^*^. Because of 
this inherently "small" system size, fluctuations in granular experimental systems are rarely 
negligible and should be considered in the analysis. 

Usually there is no general theory available to describe fluctuations in non-equilibrium 
stationary systems (NESS)^"— . A NESS is obtained in general by applying a steady driving 
force to a thermostatted system of interest: that results in a macroscopic breakdown of 
time-reversal accompanied by directed currents. The typical situation, within the realm of 
molecular fluids, is by steady shearing - giving place to a transverse momentum current - or 
by contact with two thermal reservoirs at different temperatures - leading to a steady current 
of heatr^*^. Therefore, molecular fluids are characterized by non-equilibrium conditions 
associated in general to spatially non-homogeneous configurations, i.e. currents appear in 
physical space and break isotropy or other spatial symmetries. 

An entirely different situation is present in so-called active matter, which includes driven 
granular materials, as well as "fluids" of pedestrians, flocks of animals (birds or fish), suspen- 
sions of bacteria or network of filaments (e.g. actin) in the cell, etc.—"— Common properties 
of those systems are the presence of non- conservative interactions, often mediated by a sol- 
vent, and of an external energy injection which can be - in principle - spatially homogeneous 
(for instance, animal or bacteria express their self-propulsion in the bulk of the system and 
not only at its boundaries). The situation with granular materials depends on the setup: 
experiments exist where most of the grains are in contact with the external energy source 
irrespectively of their position^i"^, resulting in a spatially homogeneous NESS. 

The present paper is devoted to investigate the structure of the velocity field in an exper- 
imental setup, which realizes an homogeneous driving mechanism, already discussed in^^. 
In homogeneous molecular fluids, the velocity field is not expected to display any struc- 
ture, as a consequence of the separation of the Hamiltonian into a kinetic and a potential 
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part "H = JC + U, the former being furtherly separated in single-particle additive contribu- 
tions, /C = Xlj-^*' that gives place to a factorized Maxwell-Boltzmann kinetic distribution 
~ Y[^exp{—PKi), with no inter-particle velocity correlations and temperature On the 
contrary, granular fluids are characterized by non- conservative interactions (kinetic energy 
is dissipated in collisions), often resulting in inter-particle velocity correlations, as well as 
velocity-position cross-correlations: such effects are best appreciated in not too dilute se- 
tups, where they give origin to visible non-equilibrium peculiarities, such as for instance the 
violation of the Einstein fluctuation dissipation theorem (FDT)^. Velocity correlations in a 
similar granular setup have been previously studied in^^, measuring them in direct space and 
- therefore - focusing on small spatial distances. In this paper we shall focus on the study 
of structure factors (in /c-space), since it highlights large space-time scales and is therefore 
more suitable to compare the experimental results with theory. Nevertheless the results 
o^ are compatible with those presented here (and in our previous publication^^), in the 
observation of exponentially space-dependent velocity correlations for the case of a rough 
vibrating plate. 



Our experiments are also interesting to discriminate between two models of "granular 
thermostats" which differ for the presence/absence of a viscous drag force (per unit of mass) 
— 7bVj proportional to the velocity v,, of the i-th particle: the presence of such a term 
allows one to deflne a typical interaction time l/'~fb with the thermostat and a flnite "bath" 
temperature (see the discussion in Section riVPp . The experiment with a rough plate is 
well described by a model with 7?, > 0, and an effective bath temperature can indeed be 
measured by looking at scales larger than a well-deflned correlation length. 



This paper has the following organization: in Section |TT] the experimental setup and 
the protocols of measurements are explained. In Section IIIII the probability distributions 
of velocities are presented. The measurements of density, transverse and longitudinal ve- 
locity structure factors are discussed in Section IIVI From a comparison of those structure 
factors with a phenomenological fluctuating hydrodynamic theory, we get access to trans- 
port coefficients, bath temperature and some typical non-equilibrium correlation lengths. 
The dependence of those quantities with external parameters are investigated in Section IVl 
Concluding remarks appear in Section |VT1 
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II. EXPERIMENTAL DETAILS 




FIG. 1. Left: a logical scheme of the experiment. Center and right: two pictures of the experimental 
setup. 

In the left frame of Fig. [1] a rationale of the experiment is presented. The center and 
right frames reproduce real photographs of the setup. See^^ and^i for movies representing 
the time-evolution of the system with different experimental parameters. The energy source 
is an electrodynamic shaker, model V450 by LDS Test & Measurement, placed on top of a 
support with tunable inclination by means of adjustable feet. An aluminum mounting plate 
is placed over the head of the shaker. The container of the granular media, which is fixed 
over the mounting plate, is a circular aluminum plate with diameter D = 200 mm, with 
perpendicular borders (whose height, 20 mm, is however irrelevant for the present setup). 
The final horizontal alignment is verified, directly on top of the plate carrying the granular 
sample, by means of a spirit level with a precision of 0.02 mm/m. The alignment check is 
repeated after every single run. The container has been produced to have a flat floor and 
then it has been submitted to sandblasting procedure with corundum mixture (aluminum 
oxide) whose nominal grain size is 1190 [im, also called "grl6" : actually, such a sand is highly 
polydisperse (grains as small as 200 [im can be found in it), and the rough surface created in 
the sandblasting process presents holes with sizes in the range [50 [im, 500 [im] , as checked 
by means of an optical microscope. Finally, a black coating is applied by electrochemical 
burnishing in order to enhance contrast for particle visualization. The shaker is fed, through 
an amplifier PAIOOO (LDS), with a sinusoidal signal at frequency /, which is converted in 
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a sinusoidal force applied to the payload. This appears as a sinusoidal movement of the 
plate with amplitude A, i.e. z{t) = A * sin(27r/t), where z{t) is the vertical coordinate 
of the plate at time t. Provided a linear response of the amplifier and of the shaker, a 
given maximum amplitude of the input signal together with a given mass of the payload 
determine a maximum acceleration z^ax = ^(27r/)^. For instance, with all other parameters 
kept constant, a change of frequency results in a change of amplitude of oscillation A ~ 1//^- 
Our experiments are done at frequencies in the range [150, 250] Hz and accelerations, in units 
of gravity acceleration g, in the range [5, 11]. We have verified that, within those limits, the 
system remains in a quasi-2d regime, with a negligible number of particles flying at heights 
larger than their diameter. Peak to peak acceleration is measured during the course of each 
run (i.e. with full payload, including moving particles) by means of an accelerometer, placed 
on the mounting plate, connected to a digital averaging oscilloscope in order to remove the 
spikes due to particle-plate collisions. 

The imaging apparatus consists of a fast camera, placed above the central point of the 
container, and a lighting dome: the former is an EoSens CL from Mikroron, the latter, placed 
at height much above the layer of particles (and therefore not interacting with them), is made 
of a led ring whose emitted light is reflected by a diffusing half-spherical surface, reaching 
the granular medium with a high degree of isotropy and uniformity. As seen in the right 
frame of Fig. Ej a very weak radial density gradient is present, so that particles in the central 
region are slightly less concentrated than those near the border. For the sake of verifying the 
importance of such a weak spatial inhomogeneity, we have performed measurements both 
in the whole container (e.g. right frame of Fig. [2]) and in a restricted region focusing on a 
central rectangular "region of interest" (roi), of size x Ly = 100 mm x 100 mm (e.g. left 
frame of Fig. |2]). Changing the optics and the distance of the camera, we manage to get, 
in both digital image with dimensions 840 x 840 pixels. Results obtained with the 

two different "roi", shown in the following sections, do not display qualitative differences 
but only weak quantitative discrepancies due to small non-homogeneities. Nevertheless, a 
larger "roi" allows one to probe modes with larger wavelengths, which is relevant in the 
comparison with theoretical predictions for structure factors at small wave-vectors. A very 
weak radial density gradient is visible in Fig. |21 due to the presence of hard boundaries 
where dissipation is stronger and density may be slightly enhanced. 

Particles used in this experiment are restricted to a single material and a few sizes, varying 
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packing fraction, i.e. varying their number A^, and shaking parameters. Here, we use spheres 
with diameter a in the range [1,4] mm made of "316-steel", which is non-magnetic in order 
to avoid interference with the shaker electromagnet. Before each experimental run, particles 
and container are washed in a ultrasound machine to eliminate those impurities which spoil 
the reproducibility of the experiment: indeed, the effectiveness of energy transfer between 
vertical vibration and horizontal spheres motion is mediated by frictional properties of the 
sphere-surface contact, which could be changed by cleanness of the apparatus. 

Steel reflectance properties are such that each particle appears as a bright and well defined 
spot, smaller than its diameter. A simple lookup table is sufficient to segment the image 
distinguishing the interesting white blobs (particles) from the black background. Particles 
appear as almost-circular blobs with an average diameter duob ^ reported in Table [T]- ranging 
from ~ 2 to ~ 8 pixels. The frame-rate facq used for image acquisition is in between [60, 250] 
frames per second. Reasons for the choice of space and time resolution are discussed below. 




FIG. 2. Two pictures from the acquisition camera (enhanced online). Left: zoom on the central 
region of interest. Right: whole system (refer also to Table [I]). The dashed box marks the zoomed 
region. 

As mentioned above, the system of interest, here, is a projection of the granular sample 
on the horizontal plane: only horizontal positions and velocities are measured. It is directly 
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observed that particles never go on top of each other, allowing to consider this system as a 
quasi-2-dimensional fluid, also in the absence of a top lid. 

The parameters used in the present paper are resumed in Table [B To have an idea of the 
relevance of interactions, we also report the mean free path A = \j{2^g2na\ where n is 
number density on the plane and g2 the value of the equilibrium pair correlation function at 
contact. A label is assigned to each experimental configuration in order to be easily referred 
to in the rest of the discussion. The rationale of the set of experiment is the following: 
experiments A-H are done with diameter 2 mm, I-T are with diameter 1 mm, U-X are with 
diameter 4 mm; experiments A-H and 1-P focus on a central region, while Q-X focus on 
the whole system; experiments A-D, 1-L, Q-T and U-X are groups of four configurations 
where only the packing fraction is changed, from 10% to 40%; experiments E-H and M-P 
are performed at constant packing fraction, but varying the shaking parameters (frequency 
or maximum acceleration). 

As a matter of fact, for all choices of parameters considered here particles perform an 
apparently random, Brownian-like motion on the plane, with low flights interrupted by 
bounces on the plate or with other particles. The study of a single particle on the plate, 
unfortunately, gives information quite different from the case of a many particle situation 
(e.g. at packing fraction > 10%), which is our system of interest. If particle's flights are not 
interrupted by collisions with other particles, repeated bounces with the plate may build up 
very high velocities (or resonances, in the case of a smooth surface)^: such a more general 
situation illustrates the fact that the particle-plate friction properties are influenced by the 
velocity of the particle itself (e.g. multiplicative noise). Nevertheless, as put in evidence 
by the comparison of our experiment with a model based on a non-multiplicative noise 
thermostat, in the many-particles configuration there is not the possibility of reaching such 
high velocities and this velocity-dependent noise regime does not take place. 

A. Acquisition of positions and uncertainty 

From each blob of pixels, the x position of the center of mass of the particle is obtained 
by averaging the average position of each row of pixels (and the same is done for columns, 
obtaining the y position). Uncertainty in the evaluation of the center of mass position is 
estimated by a Monte-Carlo procedure where discs of diameter d^oh are placed on a grid 
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TABLE L Parameters of all the experiments discussed in the present paper. 
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of unitary pixels, choosing randomly on a much finer scale (double precision, i.e. 10~^) 
the coordinate of the disc center, and then their position is computed using only covered 
pixels, with the procedure discussed before, obtaining a distribution of errors. It is seen, as 
expected, that the standard deviation of such a distribution is 5 ~ 1/ a/ duot which ranges 
- as reported in Table [T]- in ~ [0.04 — 0.16] mm, depending on the particular setup (the 
largest error is obviously associated to the experiment done by observing the full system 
with the smallest particles). It is immediately understood that even increasing the linear 
resolution by a factor 2 would not reduce significantly this error. 

B. Assessment of ballistic and diffusive timescales 

At fixed spatial resolution, an upper limit in the choice of the framerate facq for image ac- 
quisition is immediately established, given that our aim is computing velocities of particles. 
Indeed, for a population of particles with typical velocity vq ~ a/ {v'^), it is pointless measur- 
ing displacements between frames separated by a time interval At = 1/ facq < Atmin = ^ /vo, 
because such a displacement would be spoiled by the acquisition uncertainty 5 previously 
discussed. However, too large At's would be wrong as well: each particle moves ballistically 
until its fiight is interrupted by bounces with the plate or with other particles, therefore 
at large At > Atdi/f the displacement is diffusive and it does not represent the real "in- 
stantaneous" velocity anymore. For the purpose of estimating the correct time-scale Atdiff, 
we first acquire data at high frame rate (250 frames per second) and then compute the 
mean squared displacement as a function of time interval d2{t) = {[x{t) — x(0)]^): ballistic 
{d2 ~ t^) and diffusive (^2 ~ t) regimes appear clearly, marking the time-scale Atdiff which 
separates the two, see Fig. [31 We choose our final frame rate facq for the experiment such 

that Atmin < l/facg < Atdiff ■ 

It is interesting to mention that the limit of the ballistic scale is mostly independent of the 
period of oscillation of the plate, even if it may be contrary to intuition. Our main interpre- 
tation of the motion of particles on the rough vibrating surface, discussed in Section IIVDI 
is that of a Brownian motion with viscous drag determining a characteristic relaxation time 
Tf) = l/7b, basically related to the frictional properties of the sphere-surface contact, to 
which superimpose inelastic collisions with other particles occurring with a mean free time 
Tc- Such a model predicts a ballistic regime at times smaller than min[rc,rf,]: quantitative 
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FIG. 3. Mean squared displacement d2{t) for different packing fractions (j) = 
10%(H),20%(«),30%(a),40%(t), with particles' radius a = 4 mm and for driving with / = 200 
Hz (experiments with labels U,V,W,X). For all other parameters, see Table HI Error bars fall within 
the symbols. 

assessment of Tb and Tc are well consistent with the estimate of ^tdijf obtained here from 
the mean squared displacement (see discussion in Section |V|) . 

III. DISTRIBUTIONS OF VELOCITIES 

After a tiny transient regime, our system reaches a stationary homogeneous state, which is 
the subject of our research. Equilibrium statistical mechanics predicts a Gaussian (Maxwell- 
Boltzmann) distribution for velocities in a fluid. The deviation of velocities in granular 
fluids from the Maxwell-Boltzmann statistics is widely discussed in the literature^^i""— . In 
the presence of uniform driving, however, those deviations are expected to be small and well 
reproduced by polynomial corrections superimposed to the usual Gaussian statistics. 

In Fig. m two top frames, experimental rescaled velocity distributions are shown. Small 
but appreciable deviations from the Gaussian ideal curve are visible, mainly in the tails. 
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FIG. 4. Top: particle velocity distributions (rescaled to have unit variance) for experiments with 
different packing fractions (p and particles of diameter a = 1 mm (left panels) and cr = 4 mm (right 
panels), for driving with / = 200 Hz. Bottom: center of mass velocity distributions for the same 
experiments, again rescaled to have unit variance. Dashed black lines are normalized Gaussian. 
For all other parameters, see Table HI Insets: variances of the distributions, i.e. Tg and N{V^^), as 
functions of (j). 



Sonine-polynomial fits, not shown here, give a value for 02 in the range [0.01,0.05]. The 
granular temperature Tg = {1/2N) 'Ylii^f is also shown in the insets: it decreases with the 
packing fraction (p as expected due to the increasing collision rate. 

Notice that, as stated in the formula above, in the following we measure "temperatures" 
in units of squared velocities, i.e. ignoring the mass of the particle. Such a choice is 
possible because, here, temperature has only a kinetic meaning and does not correspond to 
environmental temperature (which occurs on energy scales exceedingly smaller). 
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A. Distribution of the center of mass velocity 



In the two bottom frames of Fig. HJ one can see the distribution of one component {x 
or y) of the center of mass velocity Vcm(^) = {^/N)'^f^^\'i{t), collected over different 
configurations (5 x 10^) in each experiment. It is interesting to notice how the distribution 
is well- fitted to a Gaussian. The rescaled variance N{V^^) (where by isotropy we define 
V^^ = is also plotted in the inset, as a function of the packing fraction. Notice that 

such a rescaled variance changes with the packing fraction on a much smaller energy interval 
with respect to the variation observed for Tg. 

The effect of inter-particle collisions upon the center of mass velocity is much smaller when 
compared with the velocity of the single particle: indeed collisions conserve momentum for 
each component, so that the total momentum in a region is only affected by bounces on the 
plate, and by exchanges at the boundaries of the region as well (e.g. particles going in or 
out from the region itself). The latter contribution is as smaller as larger is the size of the 
region. 

IV. STRUCTURE FACTORS 

The main goal of our study is concerned with the characterization and the analysis of 
large-scale fluctuations. The physical quantities of interest are the correlations between the 
coarse-grained hydrodynamic fields, density n{r,t), velocity u(r, t) and temperature T(r,t), 
defined as follows: 



where rj(t) is the position of particle i at time t, and the sum runs over the par- 
ticles. In particular, we focus on large spatial scale fluctuations of the system around 
the stationary homogeneous state, where the hydrodynamic fields take the values n = n, 
T = T and u = 0. We study fluctuations in momentum space and define (5a(k, t) = 



n{r,t) 




(1) 
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{(5n(k, t), 5T(k, t), M||(k, t),u±{'k, t)}, with 5a = a — a, where the Fourier transform is 

5a{k,t) = j dr 5a{r,t)e-'^'\ (2) 
and ^^(k, t) and M||(k, t) denote the shear and longitudinal modes, respectively, namely 

u±(k) = k± ■ u(k) 

«ll(k) = A;-u(k), (3) 

k± being a unitary vector such that k± ■ k = 0. In the following, we show and discuss 
the experimental data obtained for the static and dynamic correlation functions of the 
hydrodynamic fields, defined, in the stationary state, as 

= j^{Sa{k)6b{-k)), (4) 



and 



SabiKt) = ^{Sa{k,t)Sb{-k,0)). (5) 



The study of these quantities provides important information on the structural properties 
of the system and on the characteristic lengths and times governing macroscopic fluctuations. 
In particular, at equilibrium and at small k, the static structure factors, Eq. (j4j), are related 
to thermodynamic quantities^. From the measure of dynamic structure factors, Eq. ([5]), one 
can extract information on heat and sound hydrodynamic poles, and thus on the transport 
coefficients which rule thermal diffusion and propagation of pressure waves in the system. 

From the procedural point of view, we measure structure factors Sab{k,t) by directly 
applying definitions (E])-© on a set of wave- vectors k, averaging over several (typically 
5 X 10'^) configurations; subsequently, we compute sphericized averages Sab{k,t) over shells 
with average modulus k, assuming isotropy. The smallest value of k considered is kmin = 
2tx /{D /2) accordingly to the largest linear size of the system (we recall that D is the diameter 
of the container). In the case where we focus on the central region of width = D /2, the 
minimum value of wave- vector is kmin = 2tt/{Lx/2). For the largest value of k we have 
chosen a value, typically smaller than the diameter of the particles, convenient with respect 
to the relevant information we want to highlight (e.g. for the velocity structure factor it is 
useful to give a convincing evidence of the reached "plateau" at large k). 
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A. Density structure factor 

From Eq. we get the density structure factor 

^„„(k) = l(5n(k)5n(-k)) = ^ exp(-^k- (r, - r,))^ , (6) 

where the average is taken over the stationary state. 

3p ' \ ' \ ' 1 1 3^ • I • I ' 




FIG. 5. Density structure factors at various packing fractions 4> for particles with diameter a = 4 
mm and for driving with / = 200 Hz (experiments U,V,W,X). The dashed hues represent the ana- 
lytical prediction for elastic hard disks^i. In the graph on the left the experimental results are plot- 
ted versus ka, while in the graph on the right they are shown versus kcTeff with aeff/cr = cos (vr/G). 
The cartoon in the second graph illustrates the meaning of o"e// which is the effective diameter 
seen by particles which collide at an average height ~ cr/2. For all experimental parameters, see 
Table E 

In Fig. [5]we show Snn{k) measured in our system, together with the analytical prediction 
for an elastic system of hard discs. The comparison with the elastic system shows two main 
differences. First, at small k the structure factors displays an enhancement, signalling long 
range density correlations due to inelasticity. Those correlations are not the main objective 
of our analysis here and have been addressed before, see for instance^^'^'^. Notice that 
such correlations at large wavelengths do not depend monotonously on 0: higher packing 
fractions, as a matter of fact, enhance both the aggregating role of dissipation (due to 
inelastic collisions) and the dispersive role of excluded volume. Therefore, they may have 
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negative or positive influence on the density structure. Second, at larger k the theoretical 
curves for elastic hard discs are in excellent agreement with our experiment, when the 
abscissa is rescaled by a factor cos(7r/6), as shown in the right frame of Fig. O Such 
a rescaling is equivalent to consider an effective diameter of particles a^ff ~ cos(7r/6)cr. 
Indeed, particles typically collide at different heights from the plate and therefore meet at 
a distance smaller than their diameter. Such an observation is consistent with an average 
difference of height Ah ^ a/2. 

B. Transverse velocity modes 

As discussed below, in the realm of linear hydrodynamics, transverse velocity modes are 
decoupled from the other hydrodynamic fields and their behaviour can be studied separately. 
From Eq. (jlj) one has 



In Fig. [HI S±{k) is reported for different values of </>, at fixed amplitude and frequency 
of oscillation. In all cases, we find two plateau, at large and small k. The former coincides 
with the granular temperature Tg of the system and is determined by the packing fraction. 
The latter is much less dependent on and is approximately constant within experimental 
errors. Its value is expected to be mostly related to the energy injection mechanism, namely 
to the shaking parameters. Indeed, in Figs. Hand [H] the structure factor of transverse modes 
is shown upon varying the frequency and amplitude of oscillation, respectively, keeping 
fixed the packing fraction. In these cases a strong dependence of S±{k) at small k on these 
parameters is observed. At small k we have also reported (see dashed lines) the values of 
N{V^^), whose relevance at large scales is discussed below. 

The main feature arising from the figures illustrated above is a "sigmoidal" shape of the 
structure factor that can be clearly observed in the whole range of parameters. Note that 
we have intentionally avoided a too high resolution in k, to neglect "fine" effects which have 
a less clear interpretation: for instance, it would be interesting but beyond the scope of our 
analysis here, to give an interpretation of the oscillations at large k that can be observed with 
a higher fc- resolution, see Fig. [9] for an example. The main "sigmoidal" shape is remarkably 
different from what is observed in elastic fluids, where static velocity structure factors are 
flat. The profile we observe connects two regions in k space with different energy scales: one 




(7) 
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associated with the granular temperature and the other one related to the external driving 
mechanism. As it will be thoroughly discussed below, this behavior suggests the existence 
of a characteristic length-scale in the system, marking the passage from one region to the 
other. 

In order to study the characteristic times in our system, we have also measured the 
dynamic structure factor 



that is reported in Fig. [TUl for different values of the experimental parameters. Notice that, 
at small k, the decay in time of S±{k,t) is nearly exponential. Moreover, the decorrelation 
time is greater for smaller k. A more detailed analysis of such behavior is postponed to 
Section |Vl 

In Fig. [11] the probability distribution functions (rescaled to have unitary variance) of 
the transverse velocity modes are shown. For different values of the parameters, we always 
observe Gaussian distributions. 

C. Longitudinal velocity modes 

The structure factor of longitudinal velocity modes 



is shown in Fig. [121 ioT different values of the experimental parameters (see caption). The 
same qualitative behaviour is observed as for transverse modes: two energy scales are present 
at large and small k, representing the granular temperature and the energy scale associated 
to the fluctuations of the center-of-mass velocity, respectively. Again, a characteristic length 
scale can be associated with the spatial correlation of longitudinal modes, as discussed in 
more detail below. 

The temporal decay of time correlations of longitudinal modes >S'||(k, t) = -^(u||(k, t)u||(— k, 0)), 
is shown in Fig. [131 Here we find a behaviour much more complex than a simple exponential 
decay. This is due to the coupling of the longitudinal modes with the other hydrodynamic 
fields. 




(8) 




(9) 
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FIG. 6. Transverse velocity structure factors at various (p (see the legend), with particles' diameters 
£7 = 2 mm (top left), a = 1 mm (top right and bottom left) and o" = 4 mm (bottom right), for 
regions of interest "center" (top panels) and "whole" (bottom panels), and driving with / = 200 
Hz. For all experimental parameters, see Table [H Dashed lines mark the value of N{V^^) and Tg 
at low and large values of k respectively. The continuous curves are fits with formula (fT6|) . 

Finally, in Fig. [14] the probability distribution functions of M||(k) are reported. Again, we 
find that such distributions are nearly Gaussian. 

D. Comparison with theory and simulations: "bath temperature" and correlation 
lengths 

In order to discuss a theoretical model for our system, we show, starting from the defini- 
tion of S±{k), see Eq. ([7]), what the physical meaning of the two plateau appearing in Fig. [6] 
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FIG. 7. Transverse velocity structure factors at various frequencies / and fixed packing fraction 
(j) = 30% and acceleration, with particles' diameter a = 2 mm (left panel) and a = 1 mm (right 
panel). For all other experimental parameters, see Table HI Dashed lines mark the value of Niy^^ 
and Tg at low and large values of k respectively. The continuous curves are fits with formula ()16p . 





FIG. 8. Transverse velocity structure factors at various maximum acceleration of the plate Zjyiax 
and constant frequency / = 200 Hz, at packing fraction (j) = 30% and with particles' diameter 
£7 = 2 mm (left panel) and a = 1 mm (right panel). For all other experimental parameters, see 
Table HI Dashed lines mark the value of Niy^^ and Tg at low and large values of k respectively. 
The continuous curves are fits with formula ()16p . 



18 









1.1 


\ 

T 
\ 

1 \ 

\ 


— U ((t)=10%) 
■■■■ V (it)=20%) 
-- W((t)=30%) 
-- X((t)=40%) 




1-05 

g 


\ « 

- U 

w 




1 






0.95 







lo' 



FIG. 9. Transverse velocity structure factors, rescaled by T^, at various packing fractions with 
particles' diameter cr = 4 mm and for driving with / = 200 Hz (experiments with labels I, J, 
K, and L), at high values of k and finer /c-resolution, in order to show the oscillations, which are 
similar to the ones observed in the density structure factor. For all experimental parameters, see 
Table E 

is. At large values of k in Eq. ([7]) only the diagonal term survives: 

i 

whereas at small values of k we have 

^±(^^o) = ^((^l^Vij ) = y(vL) = iv(0- (11) 

For large values of k the velocity structure factor is proportional to the average square 
velocity of particles, i.e. T^, while at small k values one has Sxilz) = A^(V^„), with Ycm 
almost independent of (momentum conserving) collisions and affected only by the interaction 
with the vibrating plate. Let us note that, since in our system the momentum is not 
conserved, we have that the value of S±{k = 0) is not constrained to a 6{k) and can be 
compared to the values of the function in the neighbourhood of A; = 0. 

A model theory that fully accounts for the dissipative interactions among hard particles 
and between particles and the vibrating vessel is beyond our scope^"— . We are interested 
into an effective model that can reproduce the behavior of collective modes at large length 
scales. 
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FIG. 10. Temporal decay of the transverse velocity modes, S±{k,t), rescaled by its value at time 
t = 0, in experiments with different packing fractions (p = 10% (I), (p = 20% (J), (j) = 30% (K) and 
<j) = 40% (L), with particles' diameter a = 1 mm and for driving with / = 200 Hz. For all other 
experimental parameters, see Table HI The values of ka are the following: 0.089, 0.14, 0.2, 0.26, 
0.32, 0.51, 1, 2, 4, 8, 16. Error bars are omitted to improve legibility. 

For instance, let us consider the shear modes, which, according to linear hydrodynamic 
theory, are decoupled from the other modes of the systen>2^. The presence of an exponential 
decay for {u±{k, t)u±{—'k, 0)), which can be seen in the panels of Fig. (TUl and the measure of 
a Gaussian distribution for the real (and imaginary) part of the variable u±{k), see Fig. [TT| 
suggest that the dynamics of u±{k) is ruled by a simple linear Langevin equation. Similarly, 
for the longitudinal modes u\\{k) a Gaussian distribution for the real (and imaginary) part 
is found, as shown in Fig. [TH while the decay in time presents an exponentially damped 
oscillation. This behaviour is consistent, as discussed below, with a linear system of Langevin 
equations which couples longitudinal modes to temperature and density modulations, as 
dictated by linear hydrodynamics. 

One of the central points to be addressed in order to give a theoretical description of our 
system is the energy injection mechanism. Indeed, in a system of particles with dissipative 
interactions, the collective behaviour strongly depends on the kind of energy injection, as 



20 




FIG. 11. Probability density functions of transverse velocity modes (real part), rescaled to have 
unit variance, for experiments with different packing fractions (p = 10% (I), (f) = 20% (J), (p = 30% 
(K) and <j) = 40% (L), with particles' diameter a = 1 mm and for driving with / = 200 Hz. For all 
other experimental parameters, see Table HI The values of ka are the following: 0.089, 0.14, 0.2, 
0.26, 0.32, 0.51, 1, 2, 4, 8, 16. 

it is clear comparing, for instance, the low k tails of velocity structure factors in Refs.— 
where different kinds of model thermostats are studied. In particular, we assume that the 
velocity of the ^-th particle Vj of our experiment evolves in time according to the stochastic 
equation: 

v.(t) = -w(t) + C.(t) + F,(t), (12) 

where the force Fj(t) accounts for the inelastic collisions with other particles and we assume 
that the interaction between the vibrating vessel and each particle can be modeled as an 
equilibrium thermostat with a viscous drag — 7bVj(t) and a white noise satisfying the 2nd 
kind FDT, namely we assume for the noise: (Ci(^)Cj(^')) = '^IbTb^it — t')6ij, with ■jb the 
friction coefficient. According to this model, the energy of the center of mass of the system 
and the low k value of S±{k) are related to the temperature of the thermostat: 

SAk)^N{V^J=n. (13) 
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FIG. 12. Longitudinal velocity structure factors at various (j) (experiments with labels I, J, K and 
L). For all experimental parameters, see Table HI Dashed lines mark the value of N(y^^ and Tg 
at low and large values of k respectively. The continuous curves are fits with formula ([20]) . 



To explain the shape of the structure factors, a theory for macroscopic fields has to 
be developed from the microscopic equation (1121) . Hydrodynamics is the proper candidate 
when separation of scales exists between "fast" relaxing microscopic modes and "slow" 
relaxating macroscopic ones, a condition always to be verified in granular fiuids^'i^, but well 
satisfied in our experiment. Fluctuating hydrodynamics is a correction to such a theory when 
macroscopic scales are not large enough to neglect fluctuations. A more detailed derivation 
of the stochastic equations for the hydrodynamic fields starting from microscopic dynamics 
can be found in^. Here, we just provide phenomenological arguments for the Langevin 
equation that governs the shear modes. According to the observation of two relevant energy 
scales in our system, we assume two sources of hydrodynamic noise of different amplitude: 
an internal one, ?7/(/c, t), due to the collisions between particles, and an external one, riE{k, t), 
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FIG. 13. Temporal decay of the longitudinal velocity modes, 5||(/c,t), rescaled by its value at time 
t = 0, in experiments with different packing fractions (p = 10% (I), (p = 20% (J), = 30% (K) and 
(j) = 40% (L), with particles' diameter a = 1 mm and for driving with / = 200 Hz. For all other 
experimental parameters, see Table HI The values of ka are the following: 0.089, 0.14, 0.2, 0.26, 
0.32, 0.51, 1, 2, 4, 8, 16. Error bars are omitted to improve legibility. 

due to the action of the external thermostat. These noises have zero average and variance: 



where p is the kinematic viscosity. Then, in agreement with a 2nd kind FDT principle, 
we introduce a frictional term proportional to each source of randomness, obtaining the 
following Langevin equation: 



Such an equation accounts for all the experimental observations on shear modes presented 
above, namely the exponential decay of the time autocorrelation, the Gaussian distribution 
of Re[u±{k,t)] (and its imaginary part), and in particular the shape of the structure factor. 



iViik, t)r]j{k', t')) = 2Tgvk'^5{t - t')5{k + k') 
{VE{k, t)r]E{k', t')) = 2-i,n5{t - t')5{k + k') 
{VE{k,t)r]i{k',t')) = 0, 



(14) 




(15) 
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FIG. 14. Probability density functions of longitudinal velocity modes (real part), rescaled to have 
unit variance, for experiments with different packing fractions (p = 10% (I), (f) = 20% (J), (p = 30% 
(K) and <j) = 40% (L), with particles' diameter a = 1 mm and for driving with / = 200 Hz. For all 
other experimental parameters, see Table HI The values of ka are the following: 0.089, 0.14, 0.2, 
0.26, 0.32, 0.51, 1, 2, 4, 8, 16. 

We can say that each mode is at equilibrium, but with a different temperature. The system 
as a whole is of course out of equilibrium as long as Tg ^ Tf,. Such an observation gives 
us the way to discuss one of the most intriguing feature of the non-equilibrium structure of 
velocity correlations, namely the finding of a correlation length not present at equilibrium. 
Indeed, the temperature spectrum of modes, which is nothing but the velocity structure 
factor, can be obtained from Eq. ( fT5l) and reads as 

S4k) ^ N-Hlu^m ^ 25|±^ ^ T, . (16) 

which defines the correlation length for transverse modes 

e± = V^b. (17) 

At equilibrium, namely when collisions are elastic and Tg = Tb, the equipartition between 
modes is perfectly fulfilled and the structure factor becomes flat, i.e. we have a single 
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temperature for all the modes: S±{k) = Tb. Differently, in our non-equilibrium system of 
interest, where Tg ^ Tf,, equipartition is broken and we find an explicit dependence of S±{k) 
on k. Performing the 2D inverse Fourier transform T~^, in real space Eq. ( I16p acquires a 
plain meaning: 

r^'S^ = Tg6^'\v) + in - Tg f'^li^'-\ (18) 

where Kq{x) is the 2nd kind modified Bessel function that, for large distances, decays ex- 
ponentially 

The existence of two different and well defined energy scales has a neat consequence both in 
real and in Fourier space. In real space the distance from equilibrium — Tg determines the 
amplitude of non equilibrium correlations in the velocity field, see Eq. (1181) . In Fourier space 
a division in two subsystems that respectively equilibrated at high or low temperature takes 
place. By looking at Fig. [6] we may identify the characteristic momentum k* corresponding 
to the characteristic length ^_|_ as the conventional point which separates the subset of "hot" 
modes from the subset of "cold" modes. 

For the interested reader, a comprehensive discussion on the relations between noises, 
covariances and large scale correlations for the model ( fT2|) is addressed in^. Clearly, the 
linearization of hydrodynamic equations is only possible in a certain range of densities, and, 
in particular, when the granular fluid becomes too dense, new couplings appear and non- 
linear effects become relevant. For instance, at packing fractions where a glassy-like arrest 
of dynamics occurs^i^, the linearized hydrodynamics is no more expected to be effective. 
For the longitudinal velocity modes, we find an exponentially damped oscillating function. 
This is consistent with the fact that even within linear hydrodynamics they are coupled 
to temperature and density fluctuations and hence their temporal decay has more than 
one characteristic time-scale. Their frequency power spectrum S\\{k,u) has been calculated 
analytically for the same energy injection mechanism of Eq. (fT2!) in^. Relying on the 
fluctuating hydrodynamic equations, the static longitudinal velocity correlator turns out to 
be well approximated at low k values by the functional form: 

s\\{k) = n- in - Tg)iik^ + o{k') = Tg + + o{k'), (20) 

where does not depend only on the drag coefficient 7^ of the thermal bath and on the 
longitudinal viscosity ui but also on other microscopic parameters, as reported explicitly 
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in^. As can be seen from Fig. [121 the sigmoidal shape of S\\{k) obtained from Eq. (120|) 
is in good agreement with experiments. The fact that the low momentum approximation 
of S\\{k) recovers also the large k asymptote is not surprising: S\\{k) represents the power 
spectrum of a projection of the velocity field and is therefore strictly bounded by the energy 
scales dominant al large and small k. 



V. TIME-SCALES, TEMPERATURES AND CORRELATION LENGTHS 
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FIG. 15. Slopes of exponential decays at different values of k for S±{k,t) in experiments with 
different packing fractions = 10% (V), (p = 20% (V), = 30% (W) and (j) = 40% (X), with 
particles' diameter cr = 4 mm and for driving with / = 200 Hz. For all other experimental 
parameters, see Table HI 

From Eq. ( ITSl) the theoretical prediction for the time autocorrelation is 

S^{k,t)r^exp{-t/T{k)) (21) 
T{k)-^ = 76 + ue. (22) 

The decay in time of S±{k, t), as reported for instance in Fig. [TUl can be fit by exponentials, 
accepting a relative error in the order of 5 ^ 10% after a least squares procedure. For a 
few experiments we have analyzed the results of such fits, finding good agreement with the 
parabolic behaviour, see Fig. [151 From the parabolic fit we extract both the kinematic 
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FIG. 16. Correlation lengths rescaled by the mean free path A as function of the packing 
fraction 0, for the transverse velocity fields for different groups of experiments, with labels Q, R, 
S and T (diameter a = 1 mm) and with labels U, V, W and X (diameter cr = 4 mm), for driving 
with / = 200 Hz. The method of computing is discussed in Sec. |Vl For all other experimental 
parameters, see Table [H 
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i^_l_(mm^/s) 


(mm) 


10% 


1.1 ±0.1 


170 ± 20 


11.8 ±2 


20% 


1.1 ±0.1 


145 ± 20 


11.5 ±2 


30% 


0.97 ±0.1 


156 ± 20 


12.5 ±2 


40% 


0.96 ±0.1 


145 ± 20 


12.3 ±2 



TABLE II. 7fe, i'± and ^± as measured from exponential decays, formula (j22p . in experiments with 
particles' diameter cr = 4 mm, for driving with / = 200 Hz and all the other parameters are fixed 
to the values corresponding to label U of Table HI In particular, the dependence of ^^/A on cp is 
shown in Fig. [TBI 

viscosity u and the drag coefficient of the vibrating vessel 7b and obtain the length-scale 
at various packing fraction. Those values, for some of the experiments, are reported in 
Table [m 

Both 7b and z/^ remain roughly constant when the packing fraction is varied from 10% to 
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40%. A constant value for 7^ is in fair agreement with the interpretation that it represents 
an effective viscous drag of the particle-plate interaction, as it appears in the microscopic 
"thermostatted" model, Eq. (IT^ . The kinematic viscosity, on the other side, could weakly 
depend on 0, as well as on the internal temperature Tg, but variations of those parameters 
are not large enough to result in appreciable variations of z/_|_. Estimates, expected to be 
reliable in the dilute limit, give z/_l ~ T^Tc where Tc is the mean free time between collisions: 
the comparison with the values of Tg (reported below) is compatible with a mean free time 
Tc ~ 0.1 4- 0.5 seconds, which, in turn, agrees with the typical crossover time from ballistic 
to diffusive behavior, see Fig. [31 It is also important to notice that the order of magnitude of 
time-scale separation between the single particle relaxation 7^ and the collective (collisional) 
relaxation Tc is T^.'-fb ~ 0.1 ^0.5: this is reflected in a not very large difference between (or 
its estimate N{V^^)) and Tg. Such lack of good separation of scales makes difficult to use 
linearized hydrodynamics for quantitative predictions of velocity correlations: its predictive 
power is fully appreciated only when time-scales are well separated, as it happens with 
numerical simulations^^i^. Here, the fits of structure factor data realized with formula ( fT6ll 
and shown in Figures [61I71I8] and [T2l lead to estimates of the physical parameters not always 
consistent with those obtained independently (e.g. in the procedure discussed above for 7^ 
and i'±). For this reason in the following we have used more rough estimates of the main 
observables when comparing different experiments, as it is done in Table lllli 

The values for the correlation lengths are simply obtained by interpolating the wavelength 
where the structure factors cross half of their way between their maximum and minimum 
value, i.e. = 2Ti/k* and S±{k*) ^ s^{k^o)+s^{k^oo) ^ analogously for ^\\. Such a 
definition does not coincide in absolute value with that, coming from the theoretical model, 
■C-L ~ \/^/lb '■ indeed the latter is smaller of about a factor 3 (compare Table [III with 
Table [TTTl last four lines). Nevertheless, when comparing the same protocol of measurement 
in different experimental conditions, we expect the choice of protocol to be not too relevant. 

The granular temperatures Tg are compatible with thermal velocities which are 

about 1/5 ^ 1/10 of the maximum vertical velocity of the plate, independently of the mass 
of the particles. Fluctuations of center of mass velocity, A^(V^^), remain more stable than Tg 
when the packing fraction is changed, while they do show clear variations when the shaking 
parameters are modified: this behavior is in good agreement with the interpretation of such 
a quantity as the estimate of the "bath temperature" Ti,. Apparently, changing the vibration 
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371 
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16 
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368 


704 
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40 
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324 


691 
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30% 


796 


959 
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16 


17 


39 
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160 


370 


15 


17 


37 
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492 


892 


13 


17 


31 


40 
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355 


15 


14 


36 
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10% 
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13 
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10 
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421 
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14 


12 


17 


16 
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30% 


389 


701 


14 


21 


33 
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316 


711 
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23 


U 


10% 


536 


666 


22 
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381 


653 
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11 


w 


30% 


318 


629 


29 
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17 


20 


X 


40% 


282 


620 


32 


32 


33 
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TABLE III. Temperatures and correlation lengths measured in the experiments, following the 
phenomenological interpretation given in Section [IV D[ Lengths are measured in millimeters, time 
is measured in seconds. 
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frequency (at constant acceleration) has a more deep effect on Tb with respect to a variation 
of the acceleration at constant frequency (compare experiments E-F vs G-H, using C as a 
reference, as well as experiments M-N vs 0-P, using K as a reference). 

A general comment concerns the fact that the best agreement with the phenomenological 
fluctuating hydrodynamic model is achieved for the experiments with the fully photographed 
system (i.e. without excluding the border regions), experiments from Q to X. In those cases, 
indeed, one can reach smaller values of k both in comparison with the mean free path A 
and with the correlation lengths. Moreover, in this case the estimate of Th through N{V2n) 
is more reliable because the system under observation is really closed, i.e. center of mass 
fluctuations are not influenced by collisions of particles inside the system with particles 
outside the system. 

The correlation lengths separating the "cold" (inelastic) scales from the "hot" (elastic) 
ones, do not change appreciably as the packing fraction is changed and all other parameters 
are fixed, but they display a pronounced increase with the packing fraction when normal- 
ized by the mean free path, as can be clearly seen in Fig. [161 The meaning of such an 
observation is that the number of collisions between correlated particles increases with the 
occupied volume. The longitudinal correlation length ,^|| is - on average - slightly larger 
than the transveral one as expected also from a more detailed analysis of fluctuating 
hydrodynamics^. 

VI. CONCLUSIONS 

We have reported experiments on the measurements of density and velocity structure 
factors in a homogeneously driven quasi-2D granular material. We compared the results 
with a recent fluctuating hydrodynamic theory derived from a thermostatted microscopic 
model^. To the best of our knowledge, together with a previous publication which included 
some of the results presented here^^, this is the first systematic comparison between a gran- 
ular experiment and granular fluctuating hydrodynamics. Considering the approximations 
involved in the theory, it is remarkable to observe agreements in so many aspects: 1) the 
basic shape of the velocity structure factors, which interpolates between a "high energy" 
plateau at small k and a "low energy" plateau at large k; 2) the decay of autocorrelation of 
velocity modes, which is roughly exponential for transverse modes and damped oscillating 
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for the longitudinal ones; 3) the Gaussian character of all fluctuations of modes; 4) the 
closeness of the high energy plateau to the variance of the fluctuations of the center of mass 
velocity; 5) the fact that such a plateau strongly depends on the shaking parameters and 
much more weakly on the packing fraction. 

There are two basic assumptions underlying the theory, which are indirectly confirmed 
by the good agreement between the experiment and the predictions of the theory: 1) a 
microscopic model which includes a "thermostat" with temperature Ti, and viscous drag 
7b; 2) a sufficient separation of scales such that external noise (from the thermostat) and 
internal one (from collisions) are separated and the latter is well reproduced by a "local 
equilibrium" assumptions^. The "bath temperature" Tf, , indeed, seems to be well measured 
in the experiment by evaluating the high energy plateau of the velocity structure factors 
(at low k) or the variance of the fluctuations of the center of mass velocity multiplied by 
A^. However, the separation of scales is not significant enough to get a perfect comparison 
for the complete structure factors (e.g. with formula (fT6l) or (120|) ). The theory allows us to 
define and measure a non-equilibrium correlation length as well, which increases (in terms of 
particle mean free paths) when the packing fraction is increased, representing the separation 
of ordered "cold" scales and more disordered "hot" scales which are in equilibrium with the 
bath. 

Future investigations including a more detailed study of the vertical-to-horizontal energy 
transfer mechanism are in progress, with the aim of achieving a predictive theory for the 
properties of the effective bath Tb and 76 as functions of the shaking parameters and the 
frictional properties of the surface. 
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